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SECTION  I 


INTRODUCTION 

This  report  presents  the  results  for  the  second  year's  work  on  the 
study  of  the  dynamic  response  of  idealized  laminated  graphite-epoxy  fan 
blades  subjected  to  impact  loads.  Due  to  the  complexity  of  the  geometry 
of  the  blade,  detailed  fiber-by-fiber,  or  ply-by-ply,  "micro-mechanic" 
stress  analysis  is  not  feasible  at  this  time;  laminated  plate  theory  must 
be  used.  The  primary  purpose  of  this  study  is  to  demonstrate  the  accuracy 
of  the  equations  derived  from  lamination  theory  by  comparing  experimental 
results  with  those  calculated  from  these  equations. 

Both  in-plane  and  shear-bending  impact  results  are  included  in  the 
Report.  The  experimental  phase  of  the  present  study  consists  of  impacting 
the  edge  of  the  specimen  plate  with  a striker  plate.  The  experiment  is 
designed  to  permit  measurements  of  a one-dimensional  response.  The  velocity 
of  the  impact  is  maintained  at  a low  level  such  that  the  response  is  purely 
elastic.  The  pulse  shape  of  the  impact  loading  function  is  approximately 
rectangular.  The  duration  of  the  resulting  impact  is  controlled  by  varying 
the  striker  plate  length.  The  strain  histories  are  used  for  comparison  with 
the  analytically  predicted  wave  velocities  and  strains.  Additional  gages 
are  mounted  on  the  plate  specimens  to  monitor  simultaneity  of  impact  and 
uniformity  of  strains  across  the  plate  thickness  as  the  wave  propagates  into 
the  plate.  Three  specimen  plates  are  tested  and  the  results  compared  to 
analytical  results.  One  specimen  plate  is  fabricated  from  aluminum,  and  the 
other  two  from  graphite-epoxy,  one  is  layed-up  as  a cross-ply  and  the  other 
as  an  angle-ply.  The  aluminum  plate  has  been  used  to  ascertain  the  accuracy 
of  the  experimental  set-up  since  wave  propagation  in  aluminum  structures  is 
reasonably  well  understood. 
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The  laminated  plate  equations  used  in  this  report  are  those  derived  by 
Whitney  and  Pagano  [2].  As  discussed  in  [3],  this  set  of  equations  is  more 
accurate  than  that  of  [4],  and  may  be  considered  as  a representative  "first 
order"  laminated  plate  theory.  Several  analytical  and  numerical  solutions  to 
the  laminated  plate  equations  have  been  presented  recently.  Wang  and  Tuckmantel 
[5],  and  Moon  [6]  have  discussed  the  wave  surfaces  based  on  the  laminate!  plate 
equations  due  to  abrupt  change  in  stress.  Chow  [7]  applied  the  Laplace  trans- 
form technique  and  calculated  the  deflection  of  a laminated  plate  under  a con- 
centrated load.  Moon  [8]  solved  the  problem  of  one-dimensional  wave  propaga- 
tion due  to  a prescribed  line  load  by  the  fast  Fourier  transform  technique. 
Recently,  Sun  and  Lai  [9]  analyzed  the  response  of  a unidirectional  fiber-re- 
inforced layer  subjected  to  lateral  loading  by  applying  transform  techniques 
to  the  laminated  plate  equations  and  exact  equations.  Several  experimental 
studies  on  the  dynamic  response  of  composites  have  been  reported  in  Ref.  [10-13]. 
However,  the  major  emphasis  of  these  studies  has  been  to  determine  the  elastic 
moduli  and  dispersive  nature  of  composite  materials;  no  work  involving  the  tran- 
sient strain  histories  of  realistic  laminated  plates  subjected  to  impact  has 
been  presented. 

Detailed  description  of  the  general  set-up  of  this  program  can  be  found 
in  the  first  annual  report,  Ref.  14  and  will  not  be  repeated  here.  This  in- 
cludes the  governing  equations  used,  the  numerical  procedures,  the  specimen 
preparation,  the  in-plane  impact  experiment  arrangement,  and  the  results  of 
the  in-plate  aluminum  plate  impact.  Many  of  the  results  reported  here  have 
also  appeared  in  previous  quarterly  progress  reports.  Ref.  [15-17]. 

Results  of  this  study  demonstrate  that  the  laminated  plate  theory  used 
here  adequately  predicts  the  wave  velocities  and  transient  strains  in  symmetric 
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cross-ply  and  angle-ply  composite  plates  under  either  in-plane  impact  or 


transverse  shear-bending  impact. 

Also  included  in  this  Report  is  a discussion  of  our  progress  in  the  use 
of  HEMP,  a finite-difference  computer  code,  to  predict  material  responses, 
and  NASTRAN  to  predict  structural  responses.  In  addition,  we  present  pre- 
liminary results  of  an  impact  parametric  study  utilizing  NASTRAN  for  the  cal- 
culations. The  goal  of  this  study  is  to  identify  which  of  the  loading  function 
parameters  (pulse  duration,  pulse  shape,  object  kinetic  energy,  and  impulse) 
are  most  important  for  the  structural  response  of  a laminated  beam  subjected 
to  impact  loading. 


SECTION  II 


PLATE  IMPACT 

1.  IN-PLANE  IMPACT  STUDY 

The  description  of  the  experimental  set-up  and  the  theoretical  calculation 
of  the  in-plane  impacts  have  been  reported  previously  in  Ref.  14.  Here  we 
shall  present  directly  new  results  obtained  for  the  composite  specimens, 
a.  Results  for  Cross-Ply  Laminate 

Typical  results  from  the  strain  histories  monitored  on  the  graphite- 
epoxy  cross-ply  laminated  plate  are  shown  in  Figures  1 and  2.  Figure  1 
includes  the  results  for  the  76.2  mm  striker  and  Figure  2 the  25.4  mm  striker. 
The  strain  gage  numbers  are  shown  to  the  right  of  the  photographs.  Observa- 
tion of  the  strain  histories  from  gages  2,  4,  and  5 show  that  good  simultaneity 
of  impact  has  been  achieved  for  both  strikers.  The  maximum  difference  between 
wave  arrival  times  at  the  three  gages  was  approximately  6 ysec;  the  magnitudes 
of  the  strain  pulse  at  each  of  these  three  gages  were  nearly  identical. 

Gages  2 and  2R  demonstrate  that  little  bending  was  induced  by  either  of  the 
impacts.  The  average  wave  velocities  (7.71  mm/ysec  for  the  76.2  mm  striker 
and  7.70  mm/ysec  for  the  25.4  mm  striker)  measured  from  the  pulse  initiation 

at  gages  1,  2,  and  3 agrees  with  the  theoretical  velocity,  (A^/P)^  = 

2 

7.55  mm/ysec,  with  A^  = 0.713  MN/mm  and  P = 0.0125  gm/mm  for  the  laminated 
specimen.  These  three  gages  also  show  that  the  initial  pulse  has  little 
attenuation  or  dispersion. 

Comparisons  of  the  experimental  and  analytically  predicted  strain 
histories  are  shown  in  Figures  3 and  4.  The  analytical  results  are  based  on 
the  solution  of  equation  (4a)  of  Ref.  14  subjected  to  zero  initial  conditions 
and  a boundary  condition  at  the  impacted  edge  of 
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20  psec/div.  HORIZONTAL 
J87  Miiini/mni/cl  i v.  VERTICAL 
STRIKER  VELOCITY  = 2.423  m/sec. 


Figure  1.  Typical  Oscilloscope  Traces  lor  Graphite 
Epoxy  Cross-Ply  Plate  (76.2  mm  Striker). 
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20  usec/div.  HORIZONTAL 
187  umm/mm/div.  VERTICAL 
STRIKER  VELOCITY  = 2.617  m/sec. 


Figure  2. 


Typical  Oscilloscope  Traces  for  Graphite 
Epoxy  Cross-Ply  Plate  (25. A mm  Striker). 
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(1) 


eQ(0,t)  = 


2(An/P) 


r-  for  0 < t < t„ 
% — — 0 


= 0 


for  t0  < t 


where  tg  = 20.19  psec  and  Vg  = 2.423  m/sec  for  the  76.2  mm  striker  and 


tg  = 6.73  psec  and  = 2.617  m/sec  for  the  25.4  mm  striker. 


Observation  of  Figures  3 and  4 show  that  the  strain  magnitudes  of  the 

experimental  and  analytical  results  are  in  good  agreement.  For  the  76.2  mm 

striker  the  average  experimental  pulse  duration  was  35  psec  as  compared  to 

the  theoretical  value  of  20.19  psec;  for  the  25.4  mm  striker  the  experimental 

value  was  18  psec  versus  a theoretical  value  of  6.73  psec.  Much  of  this 

difference  is  due  to  the  delay  in  the  strain  gage  circuit  response.  The 

2 

experimental  shape  is  approximately  that  of  a sine  , whereas,  the  theoretical 
pulse  shape  is  rectangular.  This  discrepancy  in  shape  is  to  be  expected 
since  the  strain  gages  can  not  measure  discontinuities  in  loading  or  unloading. 


b.  Results  for  Angle-Ply  Laminate 

Typical  results  from  the  strain  histories  monitored  on  the  graphite- 
epoxy  angle-ply  laminated  plate  are  shown  in  Figures  5 and  6 for  the  76.2  mm 
striker  and  25.4  mm  striker,  respectively.  Observation  of  these  strain 
histories  from  gages  2,  4,  and  5 show  that  good  simultaneity  of  impact  has 
been  achieved  for  both  strikers.  Also,  gages  2 and  2R  demonstrate  that  little 
bending  was  induced  by  either  of  the  impacts.  Gages  1,  2,  and  3 demonstrate 
that  the  average  measured  wave  velocities,  7.54  mm/psec  for  the  76.2  mm  striker 
and  7.32  x 10^  mm/psec  for  the  25.4  mm  striker,  agree  well  with  the  theoretical 
velocity  of  7.39  mm/psec  where  A^  = 0.584  MN/mm  and  P = 0.0107  gm/mm".  These 
three  gages  also  show  that  the  initial  pulse  has  little  attenuation  or 
dispersion. 
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Typical  Oscilloscope  Traces  for  Crapliite  Epoxv 
An^le-Ply  Plate  (76.2  mm  Striker). 
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Typical  Oscilloscope  Traces  for  Graphite  Epoxv 
Angle-Ply  Plate  (25. A mm  Striker) 
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Figure  6. 


Comparisons  of  the  experimental  and  analytically  predicted  strain 
histories  are  shown  in  Figures  7 and  8.  The  analytical  results  are  based  on 
the  solution  of  equation  (3a)  of  Ref.  14  subjected  to  zero  initial  conditions 


and  the  boundary  condition  given  in  equation  (1),  with  t = 20.62  usee  and 
V = 2.55  m/sec  for  the  76.2  mm  striker  and  t^  = 6.87  usee  and  V = 2.55  m/sec 
for  the  25.4  mm  striker. 

Observation  of  Figures  7 and  8 show  that  the  magnitude  of  the  experimental 
and  analytical  strain  pulse  agree  well,  while  the  pulse  duration  and  shape 
do  not.  The  discussion  of  the  possible  reasons  for  these  discrepancies  has 
been  given  in  the  previous  Section. 

2 . Shear-Bending  Impact  Study 
a.  Pendulum  Impact  System 

The  basic  pendulum  impact  system  previously  designed  for  in-plane 
impacts  is  utilized  for  the  shear-bending  impact  system  (see  Figure  9 ). 

The  only  difference  is  that  the  specimen  plate  is  oriented  such  that  the 
striker  plate  strikes  the  specimen  along  its  top  surface  (near  the  edge) 
rather  than  along  its  edge. 

The  specimen  plate,  gage  locations,  velocity  measuring  system,  and 
strain  gages  are  identical  to  the  in-plane  impact  system.  The  only  difference 
is  that  gages  2 and  2R  are  used  to  ensure  a uniform  bending  impact.  The 
strain  traces  from  these  two  gages  should  be  exactly  opposite  in  magnitude 
(e.g.  top  surface  in  tension,  bottom  surface  in  compression). 
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b.  Loading  Function 

The  loading  function  for  the  specimen  in  the  transverse  impact  case 
depends  on  the  wave  motion  in  the  striker  and  cannot  be  specified  explicitly 
beforehand.  Instead,  the  wave  motions  in  the  striker  and  the  specimen  are 
solved  simultaneously,  with  certain  matching  conditions  at  the  contact  inter- 
face. The  striker  is  assumed  to  have  pure  in-plane  motion,  governed  by  the 
simple  wave  equation.  The  motion  of  the  entire  specimen  plate  is  assumed  to 
be  in  the  shear-bending  mode  where  the  three-dimensional  compression  waves 
in  the  region  close  to  the  impact  surface  have  been  neglected. 

The  three  second  order  differential  equations  governing  the  shear-bending 

motions  of  an  angle-ply  laminate,  Eqs.  (3c)  of  Ref.  14  require  three  boundary 

conditions  at  the  impacted  edge,  x=0.  The  simple  wave  equation  of  the  striker 

needs  one  boundary  condition  at  z=0,  where  the  wave  motion  in  the  plate  is 

along  x-direction  and  in  the  striker  along  negative  z-direction,  as  shown 

in  Fig.  10  . The  four  required  boundary  conditions  are  supplied  at  the  impact 

surface  x=0,  z=0,  by  (a)  matching  the  displacements  of  the  striker  and  the 

plate,  (b)  equating  the  total  axial  force  in  the  striker  to  the  total  shear 

force  in  the  plate,  (c)  assuming  the  moment  for  the  plate  is  equal  to  the  force 

in  the  striker  multiplied  by  half  the  striker  thickness,  and  (d)  assuming 

the  moment  M is  zero.  If  o(z,t)  and  u(z,t)  are  the  stress  and  displacement 
xy 

in  the  striker,  then,  these  conditions  are. 


u(0, t)  = w(0, t) 

o(0, t)  = Qx(0,t)/b 

Mx(0,t)  = a (0 , t ) b2/2 

M (0 , t ) = 0 
xy 

where  Q ,M  , and  M are  the  shear  force  and  moments  in  the  plate, 
x x xy 


(2) 

(3) 

(4) 

(5) 
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For  t > 0,  a general  relation  between  o(0,t)  and  u,fc(0,t)  can  be  derived. 
This  is  done  by  integrating  the  simple  wave  equation  along  characteristics 
in  the  striker  and  accounting  for  the  jumps  in  the  variables  across  the 
transmitted  and  reflected  waves  in  the  striker,  with  the  result 


a(0,  t)  = - pc 


0 


(-i)N"lv  -u  (°,t)  - Y 2(-l)nu,  (0,t  - — )' 

n-1,2..  c c0  - 


(6) 


for  liOhil  , t , 2£N 
C0  c0 


where  is  the  initial  velocity  of  the  striker  and  (N-l)  denotes  the  number 
of  completed  round  trip  traverses  of  the  wave  in  the  striker.  If  the  striker 
and  specimen  are  in  contact  during  the  first  traverse  of  a wave,  eqs.  (2)  and 
(3)  can  be  substituted  in  equation  ( 6 ) (with  N=l)  yielding 

1 

r-  Q(0,t)  = - pc  [V  -w , (0,t)],  0 < t < — (7) 

b x U U t Cq 

The  three  equations  (4.),  ( 5 ) and  ( 6 ) form  a set  of  proper  boundary 
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conditions  for  the  specimen  during  the  time  0 < t < — . Equations (3c) 

c0 

of  Ref.  14  are  no  longer  coupled  explicitly  with  the  striker  equation. 

If  the  stress  at  the  interface  becomes  zero  or  positive,  the  striker 
and  specimen  separate,  the  matching  condition  must  be  replaced  by  free  edge 
conditions,  and  the  subsequent  motion  of  the  specimen  can  be  solved  directly, 
independent  of  the  striker.  The  probable  time  for  separation  to  occur  is 
at  N2 2/Cq,  the  time  when  the  wave  in  the  striker  due  to  the  original  impact 
returns  to  the  striker-specimen  interface.  To  determine  if  separation  occurs 
at  these  instants,  the  normal  stress  must  be  found  immediately  after  the  wave 
is  reflected  and  transmitted  at  the  interface.  This  is  accomplished  by 
assuming  that  both  bodies  remain  in  contact  and  by  considering  equations  (2), 
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(3),  and  (6)  for  t = , and  the  jump  relationship  between  Qx  and 

C0 

[Q  ] = - pc-h [w  ] (8) 

X J , t 

where  the  bracket  denotes  discontinuous  values  of  variable  enclosed.  Ibis 

2NJ+  i 

is  a system  of  4 equations  at  t * — in  the  4 unknowns  o,  u>t>  w>(.,  and 

0 

0 . The  resulting  expression  for  the  normal  stress  is 
Xx 


f3© 


C0b+C3h  C 


A »<■)  t u,  C » 

3h  \ coJ  0 c0 


N-l  N 1 / , \ n L 2Ni+  2nA  I 

(-D  v + 2pc  r (-D  u,  (0,  — - — ) | 

0 ° n=l,2, . . V 0 L0  / 1 


If  the  value  of  the  right  hand  side  of  equation  (9)  is  equal  to  or 
greater  than  zero,  separation  occurs.  If  not,  the  bodies  remain  in  con- 


The  procedure  for  calculating  the  responses  of  the  specimen  when  separa- 
tion does  not  occur  after  the  first  traverse  of  the  wave  in  the  striker  is 

2i  4i 

straightforward.  Equations  (3c)  of  Ref.  14  are  solved  for  < t <:  — 

C0  c0 

subject  to  the  boundary  conditions  (4),  (5),  and  equation  (6)  with  N=2  and 

4 £ 

equations  (2)  and  (3)  inserted.  If,  after  t = — , the  striker  still  does 

C0 

not  separate,  equations  (3c)  of  Ref.  14  are  solved  in  the  time  regime 

4 £ 62, 

— < t < — subject  to  the  boundary  condition  ( 4 ),  ( 5 )>  and  equation  ( 6 ) 
c c 

o o 

with  N=3  and  equations  (2)  and  (3)  inserted.  This  procedure  is  repeated 

for  every  — in  time  until  the  striker  is  shown  to  separate  by  equation  (9  ) 
c0 

or  sufficient  calculations  have  been  achieved. 

The  procedure  for  the  cross-ply  laminate  is  similar  to  that  just  out- 


lined for  the  angle-ply  laminate. 
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c.  Results  for  Aluminum  Plate 

In  order  to  check  out  the  accuracy  of  the  shear-bending  impact 
system,  a 2024-T351  aluminum  plate,  rather  than  a laminate  specimen,  was 
tested  first.  Typical  preliminary  results  from  the  strain  histories 
monitored  on  the  aluminum  plate  are  shown  in  Figure  11.  The  gage  numbers 
are  shown  to  the  right  of  the  photograph.  Observation  of  the  strain  histories 
from  gages  2,  4,  and  5 show  that  good  simultaneity  of  impact  has  been 
achieved.  Gages  2 and  2R  demonstrate  that  we  have  achieved  uniform  bending 
across  the  thickness  of  the  plate.  The  wave  velocity,  as  measured  from  the 
pulse  initiation  of  gages  1,  2,  and  3 was  3.264  mm/psec.;  this  agrees  well 
with  the  theoretical  shear  velocity  of  3.11  mm/psec  (k  = 1.0).  The  strain 
traces  in  Figure  11  also  show  the  arrival  of  the  wave  reflecting  from  the  back 
edge  of  the  plate.  Based  on  these  preliminary  results,  we  have  concluded 
that  this  experimental  system  is  satisfactory. 

d.  Results  for  Cross-Ply  Laminate 

Application  of  eq.  (9)  to  the  impacts  involving  the  76.2  mm  striker 
showed  that  it  would  separate  after  21/ c^  while  the  25.4  mm  striker  was 
shown  not  to  separate  (Section  II-2b)  until  6£/cQ. 

Figure  12a  presents  typical  results  from  the  strain  histories  on  the 
cross-ply  laminate  produced  by  the  transverse  shear-bending  impact  of  the 
76.2  mm  striker.  The  trace  numbers  again  refer  to  the  strain  gage 
in  Fig.  1 of  Ref.  14.  At  least  four  impact  tests  with  almost  identical  strain 
history  traces  were  obtained  for  each  impact  case.  Again  the  strain  histories 
from  gages  2,  4,  and  5 show  that  good  simultaneity  of  Impact  was  achieved. 

Tlie  maximum  difference  between  wave  arrival  times  at  these  three  gages  was 
approximately  5 qsec;  except  for  gage  4,  the  magnitudes  of  the  strain  pulse 


at  these  gages  were  found  to  be  essentially  identical.  The  reason  for  the 


20 


40  usec/div.  HORIZONTAL 
400  , mm/mm/d iv.  VERTICAL 


Figure  11.  Typical  Oscilloscope  Traces  for  the  Shear-Bending 
Impact  of  an  Aluminum  Plate  (7b. 2 mm  Striker). 
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slight  discrepancy  in  the  magnitude  at  gage  4 is  that,  at  the  impact  surface, 
uniform  contact  during  the  impact  was  not  achieved,  because  of  minute  ripples 
on  the  surface  of  the  laminate.  For  instance,  for  an  impact  speed  of  2.5  m/sec, 
a 0.0254  mm  (0.001")  irregularity  on  the  impact  surface  causes  a 10  psec 
discrepancy  in  arrival  time.  Gages  2 and  2R  demonstrate  that  near  perfect 
bending  was  induced  by  the  impact. 

The  laminated  plate  properties  are  calculated  from  lamina  properties  and 
are  listed  in  Table  1.  The  striker  velocity  was  2.55  m/sec.  The  average 
wave  velocity  determined  from  the  peak  strain  at  gages  1,  2,  and  3 was 
1.67  mm/psec,  which  is  much  closer  to  the  theoretical  shear  velocity  of 
1.53  mm/psec,  than  the  theoretical  bending  velocity  of  7.69  mm/p sec,  indicat- 
ing that  the  pulse  is  primarily  a shear  wave.  There  is  some  uncertainty 
about  the  measured  shear  wave  speed.  For  instance,  if  measured  from  the  first 
disturbance  at  the  gages,  the  shear  velocity  would  be  1.81  mm/psec. 

Figure  13  shows  that  the  pulse  shape  and  magnitude  of  the  experimental 
results  are  in  general  agreement  with  those  from  the  theoretical  calculation. 

The  analytical  strain  at  25.4  mm  contains  a precursor  bending  wave  travelling 
with  the  plate  wave  velocity,  as  indicated  by  the  abrupt  wave  front.  (The 
precursor  at  127  mm  is  not  plotted).  The  strain  gage  did  not  pick  up  this 
precursor,  probably  due  to  its  limitation  in  frequency  response.  However, 
the  fact  that  the  bulk  of  the  energy  is  travelling  with  a velocity  nearly 
equal  to  the  shear  velocity  is  evident. 

The  experimental  strain  histories  of  the  cross-ply  laminate  subjected 
to  impact  by  the  25.4  mm  striker  were  similar  in  nature  to  those  of  the 
76.2  mm  striker  and  are  shown  in  Fig.  12b.  The  previous  discussion  of  the 
experimental  results  also  apply  to  this  impact  situation.  Figure  14  presents 
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Figure  14.  Theoretical  and  Experimental  Strain  Response  of  Cross-Ply 
Laminate  Due  to  Transverse  (Shear-Bending)  Impact 
(25.4  mm  Striker). 
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a comparison  of  the  25. A mm  striker  induced  strain  histories  on  the  cross-ply 


laminate  with  analysis.  We  note  general  agreement  between  the  experimental 
and  analytical  pulse  shapes  and  pulse  magnitudes.  'Hie  average  velocity  based 
on  first  disturbance  was  1.90  mm/u sec  and  based  on  peak  strain  was  1.58  mm/usec. 

The  theoretical  shear  velocity  for  the  cross-ply  was  1.53  mm/usec. 

e.  Results  for  Angle-Ply  Laminate 

Results  of  the  angle-ply  specimen  subjected  to  transverse  impact  by  the 
76.2  mm  striker  are  shown  in  Figure  15.  The  striker  impact  velocity  is 
2.55  m/sec.  Application  of  eq.  (9)  to  the  impacts  involving  the  76.2  mm 
striker  showed  that  it  would  separate  after  25/cq.  The  average  of  the 
measured  wave  velocities,  based  on  first  disturbance  at  gages  i,  2,  and  3 
was  2.11  mm/usec.  The  average  velocity  based  on  peak  strain  was  1.39  mm/. sec. 

The  theoretical  shear  velocity  for  the  angle-ply  was  computed  to  be  1.58  mm/usec. 

The  agreement  between  the  theoretical  and  experimental  shear  velocities  is  not 
as  good  as  that  in  the  in-plane  velocity,  c^.  The  shape  and  magnitude  of  the 
pulses,  however,  do  show  fairly  good  agreement. 

Figure  16  presents  the  comparison  of  the  25. A mm  striker  induced  strain 
histories  on  the  angle-ply  laminate  with  analysis.  For  this  case,  equation  (9) 
predicted  separation  after  t = A(/c^.  Again,  we  note  the  general  agreement  of  j 

the  pulse  shapes  and  magnitudes.  The  average  velocity  based  on  first  distur- 
bance was  2.00  mm/usec  and  based  on  peak  strain  was  1.78  mm/usec.  The 
theoretical  shear  velocity  for  the  angle-ply  was  1.58  mm/usec. 
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Figure  16.  Theoretical  and  Experimental  Strain  Response  of  Angle-Ply 
Laminate  Due  to  Transverse  (Shear-Bending)  Impact 
(25.4  mm  Striker) 
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3.  Summary  and  Conclusions 


We  have  obtained,  both  experimentally  and  analytically,  the  transient 
responses  of  aluminum  and  graphite-epoxy  cross-ply  and  angle-ply  laminates 
due  to  in-plane  and  shear-bending  impacts.  A comparison  of  the  wave  speeds 
is  summarized  in  Table  2 and  a comparison  of  strain  histories  is  presented 
in  Figures  3,4,7,8,13,14,15,  and  16. 

Based  on  these  comparisons  between  experimental  and  theoretical  results, 
we  conclude  that  the  laminated  plate  theory  used  here  adequately  predicts 
the  wave  velocities  and  transient  strains  in  symmetric  cross-ply  and  angle-ply 
composite  plates  under  either  in-plane  impact,  or  transverse  shear-bending 
impact.  The  wave  length  to  plate  thickness  ratio  for  the  results  reported 
here  is  approximately  8. 


III.  STATUS  OF  LOCAL  AND  STRUCTURAL  RESPONSE  CALCULATIONS 


1 . HEMP  - Local  Response 

We  have  made  progress  in  modifying  HEMP  to  accommodate  anisotropic 
materials.  We  have  incorporated  into  the  code  constitutive  equations  for 
general  triclinic  material  in  plane  strain  problems,  and  for  transversely 
isotropic  material  for  axisymmetrical  problems.  With  slight  further  modi- 
fication, the  case  of  plane  stress  can  be  handled.  We  are  in  the  process 
of  debugging  these  modifications  by  running  a sample  problem.  The  follow- 
ing is  a summary  of  the  equations  used  in  the  anisotropic  constitutive  re- 
lations . 

In  as  much  as  HEMP  traces  the  position  of  the  particle,  it  can  be 
considered  a Lagrangian  program.  However,  the  governing  equations  such  as 
the  equation  of  motion,  constitutive  relations,  etc.  are  written  with  respect 
to  coordinate  system  fixed  in  space.  In  this  sense,  HEMP  could  be  considered 
as  a hybrid  Eulerian-Lagrangian  computer  code. 

A composite  can  be  "smeared"  and  treated  as  a homogeneous  anisotropic 
material.  The  constitutive  relation  of  HEMP  has  been  modified  to  handle  such 
cases  in  the  elastic  range.  The  new  constitutive  relations  are 
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where  S..  are  the  components  of  the  deviator  of  the  stress  tensor;  e, 
ij  kx. 

are  the  components  of  the  deviator  of  the  strain  tensor;  a is  the  spherical 
component  of  stress;  Is  the  mean  normal  strain;  and  c ^ ^ , are  the  stiffness 
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constants  relating  stress  and 
previous  state  of  stress  into 


strain.  The  term  f..  is 


the  current  configuration 


for  transforming  the 
and  is  given  by 
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where  a.^  are  the  direction  cosines  and  the  superscript  "0"  signifies  the 
previous  state.  The  dot  indicates  the  total  time  derivative.  All  indices 
in  equations  10-12  run  over  1,2,3.  The  other  governing  equations  of  HEMP 
remain  essentially  unchanged. 

Equations  10-12  have  been  programmed  to  handle  axisymmetrical  and  plane 
strain  problems.  Some  slight  modifications  to  HEMP  still  have  to  be  made  in 
order  to  handle  plane  stress  problems.  The  axisymmetrical  problems  are  limited 
to  transversely  isotropic  materials;  whereas,  the  plane  strain  case  can  handle 
up  to  a triclinic  material.  However,  the  plane  strain  case  may  later  be 
modified  to  handle  only  up  to  a monoclinic  material  in  order  to  save  computa- 
tion time. 

As  a means  of  debugging,  an  isotropic  problem  consisting  of  an  iron 
cylinder  impacting  another  iron  cylinder  has  been  run.  Early  indications  are 
that  Drexel's  version  of  HEMP  containing  the  previously  mentioned  anisotropic 
elastic  constitutive  relations  gives  good  results.  The  next  step  in  the  de- 
bugging process  will  be  to  run  a problem  of  a rigid  cylinder  impacting  a 
semi-infinite  composite  of  [0,  + 45,  + 90]  lay  up.  This  composite  can  be 
"smeared"  to  give  a homogeneous  transversely  isotropic  material.  Then  the 
problem  can  be  solved  by  HEMP  as  an  axisymmetrical  case. 


2. 


NASTRAN  - Structural  Response 
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Three  structural  impact  problems  are  presently  being  calculated  by 
NASTRAN.  One  involves  the  transverse  impact  of  a rod  on  a beam,  the  second 
involves  the  transverse  impact  of  an  angle-ply  laminated  plate  (similar  to 
the  problem  discussed  in  Section  II-2)  and  the  other  involves  a late-stage 
equivalence  study  in  structural  impact  utilizing  a laminated  beam.  The  progress 
of  each  is  discussed  in  the  following  paragraphs. 

a.  Transverse  Impact  of  Rod  on  Beam-NASTRAN  Analysis 

An  important  objective  of  Drexel's  F.O.D.  research  program  is  to  obtain 
the  transient  elastic  response  of  impacted  structures.  Unfortunately,  most 
realistic  F.O.D.  impact  problems  are  difficult  to  analyze  due  to  the  compli- 
cated geometries  and  loadings  involved.  In  addition,  these  problems  are 
further  complicated  by  the  use  of  anisotropic  materials  and  laminated  con- 
struction. Traditional  analyses  which  utilize  partial  differential  equations, 
such  as  the  theory  of  elasticity  or  beam,  plate,  and  shell  theories,  are 
practically  limited  to  simple  geometries  and  boundary  conditions.  Numerical 
methods,  especially  the  finite  element  method,  seems  to  be  the  only  recourse. 

In  the  following  sections  an  application  of  this  method  to  a two-body  impact 
problem  will  be  presented. 

F.E.M.  Computer  Code  and  Impact  Problem 

To  solve  complicated  impact  problems  using  the  finite  element  method 
a computer  code  which  possesses  the  following  characterist ics  is  necessary:  a 
large  number  of  types  of  elements,  anisotropic  material  capability,  flexible 
boundary  condition  modeling  techniques,  and  more  than  one  transient  solution 
techniques.  We  believe  that  the  NASTRAN  (Ref.  18)  computer  code  possesses 
these  characteristics. 
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In  order  to  compare  our  NASTRAN  analysis  with  existing  analytical  and 
experimental  results,  we  solved  the  problem  of  the  normal  impact  of  a semi- 
infinite isotropic  rod  on  an  infinite  isotropic  beam  (Ref.  19).  In  the  re- 
maining sections,  the  details  of  the  NASTRAN  modeling  of  this  impact  problem 
and  a comparison  with  the  results  (experiment  and  analysis)  will  be  presented. 

NASTRAN  Boundary  Condition  Modeling 

Because  the  rod  is  governed  by  the  simple  wave  equation,  the  interaction 
of  the  rod  on  the  beams  can  be  approximated  by  an  appropriate  boundary  condition 
on  the  beam.  In  addition,  the  symmetry  condition  is  used.  At  the  bar-beam 
interface,  the  local  deformation  of  the  beam  is  neglected,  only  the  beam  mode 
of  deflection  is  considered.  The  force  at  the  impact  end  of  the  bar,  which  is 
assumed  to  be  equal  to  the  shear  in  the  beam,  can  then  be  calculated  as  a product 
of  the  acoustic  impedance,  cross-section  area,  and  the  velocity  (V^  + 6,  ).  In 
terms  of  the  coordinate  system  shown  in  Fig.  17,  the  boundary  conditions  of  the 
beam  are  x = 0 

2Q  = - PcAr(VQ  + 6,t)  (13) 

and 

6=0.  (14) 

In  equation  (13),  Q is  the  shear  force,  6 is  the  deflection  of  the  beam, 

and  P,  c,  and  A are  respectively  the  density,  wave  speed,  and  cross  sectional 

area  of  the  rod. 

Equations  (13)  and  (14)  can  be  recast  into  finite  element  form  by  con- 
sidering the  impact  to  occur  at  node  1 and  by  making  the  following  substitutions. 
Replace  Q | x_q  by  an  external  force  acting  on  node  1 in  the  y-dir.,  F2^'^;  also 
replace  6 by  the  translational  degree  of  freedom  in  the  y-dir.  at  node  1, 

1/2^^;  and  replace  by  the  rotational  degree  of  freedom  about  the  z-axis 

at  node  1,  The  equivalent  finite  element  form  of  equations  (13)  and 

(14)  are  then 
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and 


2F2(1)  + pcArU2,([)  = - pcArV0 

(1) 


(15) 


U6  = 0 (16) 

For  most  finite  element  computer  codes  a boundary  condition  of  the  form 
of  equation  (15)  cannot  be  used,  however  through  the  use  of  NASTRAN's  transfer 
function  capability  equation  (15)  can  be  recast  into  an  acceptable  form.  In 
essence.  Equation  (15)  is  treated  as  the  governing  equation  of  a fictitious 
element,  with  F2 ^ and  pcArV()  as  the  two  displacements  of  the  element.  Due  to 
the  U2,^^  term,  the  viscous  damping  matrix  must  be  used  although  it  has  only  one 
non-zero  term.  This  fictitious  element  adds  two  additional  rows  and  columns  to 
all  matrices. 

To  accomplish  this,  we  first  augment  the  matrix  governing  equations  by  two 
rows  and  columns.  This  is  accomplished  by  using  an  EPOINT  bulk  data  card.  Call 
these  non-structural  degrees  of  freedom  <^qq  and  e>QQ’  anc*  take  e^^  equal  to 
F2^‘>  and  e^^Q  equal  to  pcA^V^.  Second,  from  the  finite  element  equation  for 
F2^^  in  terms  of  the  discrete  degrees  of  freedom  obtain  an  equation  which  re- 
lates e^QQ  to  the  appropriate  discrete  degrees  of  freedom.  This  equation  has  the 
following  form 


M2 (1^U2,  + M2^)U2,(2)  + M6(2^U6,  ^ + K2^1^U2^1)  + K2(2'U2(2)  + 

’tt  ’tt  tt 


+ K6(2)U6(2)  - e 


(17) 


100 


0 


where  the  M's  are  the  consistent  mass  terms  and  the  K's  are  the  stiffness 
terms.  Third,  use  e|QQ  and  e>oo  to  recast  equation  (15)  into  the  following 


2e100  + ^V2’^  + e200  " °* 

Finally  define  an  external  loading  on  e.7(^  such  that 

e cA  V . 

200  r 0 


(18) 


(H) 


The  result  of  this  procedure  is  that  F2^^  is  defined  implicitly  as  a function 
of  U2,^^.  Equations  (17)  and  (18)  are  generated  in  NASTRAN  by  use  of  TF  hulk 
data  cards.  Equation  (19)  is  generated  in  two  steps.  First  a TF  hulk  data  card 
is  used  to  generate  the  left  hand  side  of  the  equation,  and  second  T1.0AD1, 
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1 


TABLE1,  and  DAREA  bulk  data  cards  are  used  to  generate  the  load  pcA^V^.  Note 

for  quiescent  initial  conditions  the  value  of  pcA  V„  must  be  fast  ramped.  The 

r 0 r 

preceding  procedure  was  utilized  to  affect  boundary  condition  (15)  in  our 

NASTRAN  analysis. 

NASTRAN  Problem  and  Comparison  with  Experiment 

The  finite  element  geometry  used  in  the  analysis  is  shown  in  Fig.  17. 

The  discretization  was  chosen  in  order  that  we  could  accurately  include  all 
modes  up  to  the  fifteenth  symmetric  mode,  and  the  overall  length  of  the  beam 
was  chosen  so  that  the  solution  at  x = 25.4  mm  and  76.2  mm,  would  not  be  in- 
fluenced by  wave  reflection  from  the  boundary  for  at  least  100  usee. 

The  stiffness  matrix  was  formed  from  beam  elements  which  have  a shear  de- 
flection correction  and  the  mass  matrix  was  formed  >y  the  consistent  approach. 
The  transient  equations  were  solved  using  NASTRAN' s direct  integration,  rigid 
format.  'Hie  time  step  chosen  at  At  = 0.5  usee,  which  is  approximately  one 
third  of  the  period  of  the  fifteenth  mode. 

In  Fig.  18  a comparison  is  presented  between  the  NASTRAN  solution  and  the 
analytical  and  experimental  results  of  Ref.  [19],  for  the  normal  impact  of  a 
rod  and  beam.  The  NASTRAN  solution  is  in  good  agreement  with  the  experimental 
results.  The  analytical  solution  of  Ref.  [19]  was  obtained  from  numerical  in- 
version of  che  Laplace  transformation  of  the  Timoshenko  beam  equations  with 
boundary  conditions  (13)  and  (14).  The  NASTRAN  solution  does  pick  up  the 
change  in  shape  of  the  response  at  x = 1"  and  x = 3",  and  it  does  reasonably 
yield  the  peak  response  at  these  two  locations.  Oscillations  in  the  NASTRAN 
solution  may  be  eliminated  by  using  a smaller  time  step. 

On  the  basis  of  the  good  agreement  between  the  NASTRAN  solution  and  the 
experimental  data  for  this  two  body  impact  problem,  we  will  apply  this  technique 
to  more  complicated  two-body  impact  problems. 
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BOUNDARY  CONDITIONS 

1:  2F2(1)  + p cA  U2,(1)  = - p cA  Vn 

r t r 0 

U6(1)  = 0 (SYMMETRY) 

31:  U2(31)  = U6(31)  = 0 (FIXED) 


INITIAL  CONDITIONS 

):  U2(i)  = U6(l)  = 0 

U2,(1)  = U6,(U  = 0 
’t  t 


i = 1,2 31 


Figure  17.  Finite  Element  Geometry  and  Auxiliary  Conditions 
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Figure  18.  NASTRAN  SOLUTION  OF  THE  IMPACT  OF  A SEMI-INFINITE 
STEEL  ROD  ON  AN  INFINITE  STEEL  BEAM 
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b.  Transverse  Impact  of  Laminated  Plate 

The  response  of  a graphite-epoxy  [ (0/45/0/-45)^/0]s  angle-ply  laminated 
plate  (609.6  mm  x 355.6  mm  x 7.11  mm)  subjected  to  a transverse  shear  impact 
by  a semi-infinite  striker  is  being  calculated  by  NASTRAN.  The  material  and 
lay-up  of  the  striker  is  identical  to  the  impacted  plate;  the  problem  is 
identical  to  that  described  in  Section  II-2e  except  that  the  striker  is  semi- 
infinite rather  than  either  25.4  mm  or  76.2  mm  long.  Our  initial  run  utilized 
NASTRAN's  plate  element  (25.4  mm  square)  with  the  direct  transient  solution 
technique.  We  have  used  the  transfer  function  capability  to  prescribe  the 
two-body  impact  boundary  conditions  and  a time  step  of  5 psec.  For  comparison 
purposes,  we  have  solved  the  same  problem  by  applying  the  MCDU-22  code  [20] 
to  equations  (3c)  of  Ref.  14  subjected  to  the  semi-infinite  striker  impact. 
Figure  19  presents  the  comparison  of  the  transient  strain  histories  at  the 
25.4  mm  and  76.2  mm  locations.  The  results  of  this  comparison  demonstrate 
that  the  NASTRAN  computation  qualitatively  predicts  the  response  (e.g.  shape) 
but  not  the  magnitude  at  the  earlier  times.  We  are  presently  recomputing  the 
NASTRAN  strain  histories  using  12.7  mm  square  elements.  The  goal  of  this  phase 
of  the  study  is  to  assess  the  adequacy  of  NASTRAN's  plate  elements,  impact 
boundary  condition  prescription,  and  solution  technique  for  structural  impact 
response  predictions. 
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Late-Stage  Equivalence  Study  in  Structural  Impact 


C . 

This  phase  of  the  program  involves  a parametric  study  of  a 
laminated  beam  model  of  a fan  blade  subjected  to  different  impact  loadings 
on  the  end  of  the  blade.  The  purpose  of  this  study  is  to  determine  the  late- 
stage  equivalence  of  structural  responses  caused  by  impacts  defined  by 
similar  or  dissimilar  parameters  (e.g.  kinetic  energy,  impulse,  pulse  shape, 
and  pulse  duration).  Figure  20  describes  the  overall  problem  and  NASTRAN 
modelling  being  utilized  in  this  study.  Figure  21  describes  the  loading  im- 
pact parameters  used  for  each  part  of  this  study.  To  date  we  have  obtained 
1000  Msec  of  transient  response  at  the  5",  10"  and  19"  locations  of  the  blade 
model  for  each  of  loading  conditions  described  in  Figure  21  . Some  of  these 
preliminary  results  are  presented  in  Figures  22  thru  25.  Figures  22  and  23 
present  the  transient  shear  and  moment  responses  for  the  loading  conditions 
having  equal  impulses.  We  see  from  these  Figures  that  the  magnitude  of  the 
responses  due  to  the  larger  magnitude  (but  smaller  pulse  duration)  input 
pulse  is  initially  much  larger,  but,  with  time  this  difference  in  response 
magnitudes  decreases.  This  trend  is  particularly  evident  at  the  19"  location. 
Figures  24  and  25  present  the  transient  shear  and  moment  responses  for  the 
loading  conditions  having  equal  kinetic  energies  (elastic  striker  - unequal 
impulses).  For  this  case  we  see  that  the  pulse  with  the  higher  magnitude 
(lower  impulse)  produces  responses  which  are  much  smaller  at  longer  times. 
These  preliminary  results  indicate  that  the  impulse  is  an  important  para- 
meter in  any  impact  boundary  condition  specification.  We  are  presently 
computing  some  of  the  responses  for  1500  psec  to  verify  the  aforementioned 
trends  and  are  pLotting  the  results  for  the  pulse  duration  and  pulse  shape 
studies . 


42 


• MATERIAL  = GRAPHITE  EPOXY 
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Figure  20.  Late-Stage  F.quivalence  in  Structural  Impact-Problem  Description 
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Figure  22.  Transient  Shear  Response  for  Loading  Conditions  Having 

Equal  I^julses 
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